Production of ultracold NH molecules by sympathetic cooling with Mg 
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We carry out calculations on Af-changing collisions of NH ("^S") molecules in magnetically trap- 
pable states using a recently calculated potential energy surface. We show that elastic collision rates 
are much faster than inelastic rates for a wide range of fields at temperatures up to 10 mK and 
that the ratio increases for lower temperatures and magnetic fields. If NH molecules can be cooled 
to temperatures approaching 10 mK and brought into contact with laser-cooled Mg then there is a 
good prospect that sympathetic cooling can be achieved. 

PACS numbers: 34.50.-s,34.50.Cx,37.10.Mn,37.10.Pq 



There is great interest in the production of samples of 
cold molecules, at temperatures below 1 K, and ultra- 
cold molecules, at temperatures below 1 mK. Ultracold 
molecules have many potential applications in areas rang- 
ing from precision measurement to quantum computing. 
They also offer new possibilities for quantum control and 
controlled ultracold chemistry. 

There have been considerable successes in producing 
ultracold molecules in laser-cooled atomic gases [l| , both 
by photoassociation and by magnetoassociation 
Ni et al. [3| have very recently succeeded in producing 
ultracold KRb in its ground rovibrational state by mag- 
netoassociation followed by stimulated Raman adiabatic 
passage (STIRAP). However, such methods are limited 
to molecules formed from atoms that can be laser-cooled, 
such as the alkali metals. A wider range of molecules 
can be cooled directly from high temperature to the mil- 
likelvin regime, using methods such as buffer-gas cooling 
@ and Stark deceleration @. Low- field-seeking states 
of these cold molecules can then be confined in electro- 
static and magnetic traps. However, at present the low- 
est temperature that can be achieved for directly cooled 
molecules in static traps is around 10 mK. The ma- 
jor challenge in this field is to find ways to cool such 
molecules to the ultracold regime. 

One of the most promising proposals for second-stage 
cooling is sympathetic cooling, in which molecules are 
brought into contact with a laser-cooled atomic gas that 
is already ultracold. The hope is that thermalisation will 
occur to produce ultracold molecules. However, atom- 
molecule potential energy surfaces are often strongly 
anisotropic and the anisotropy may drive fast inelastic 
collisions (relaxation). Such collisions may prevent sym- 
pathetic cooling, because they release kinetic energy and 
cause trap loss. A commonly stated rule of thumb is that 
elastic collisions must be at least a factor of 100 faster 
than inelastic collisions if sympathetic cooling is to suc- 
ceed. 

In previous work, we have used high-level electronic 
structure calculations to calculate interaction potentials 
for a variety of systems that are candidates for sympa- 
thetic cooling. The systems investigated include OH with 
Rb 0, [1] and NH3 and NH (^E) with alkali-metal and 
alkaline-earth atoms 0, [l^ . Most of these systems were 



found to have interaction potentials with deep wells and 
strong anisotropy, and several of them also have ion-pair 
states that are expected to cause additional inelasticity. 
It is unlikely that sympathetic cooling would work for 
molecules in low- field- seeking states in systems with high 
anisotropy. However, Mg -I- NH and Be -I- NH were found 
to be much less anisotropic and their ion-pair states are 
likely to be energetically inaccessible in low-energy colli- 
sions ■ Cold NH in its ground state can be cooled 
in a helium buffer gas and confined in a magnetic trap 
[ill . [T^ , and laser cooling of Mg to sub-Doppler temper- 
atures has recently been achieved jl^. NH can also be 
decelerated and trapped electrostatically in its excited 
electronic state [3, [iBl and there is a proposal to 
transfer the molecules to the ground state and accumu- 
late them in a magnetic trap [l^. The purpose of the 
present paper is therefore to use quantum collision cal- 
culations to investigate whether sympathetic cooling of 
magnetically trapped NH by Mg is a good prospect. 

The energy levels of NH in a magnetic field are most 
conveniently described using Hund's case (b), in which 
the molecular rotation n couples to the spin s to produce 
a total monomer angular momentum j. In zero field, 
each rotational level n is split into sub-levels labeled by 
J. In a magnetic field, each sublevel splits further into 
2j -I- 1 levels labeled by mj, the projection of j onto the 
axis defined by the field. For the n — Q levels that are 
of most interest for cold molecule studies, there is only 
a single zero-field level with j = 1 that splits into three 



components with rrij 



-1, and —1. Molecules in the 



rrij = -|-1 state are low-field-seeking and can be confined 
in a magnetic trap, whereas those in the mj = and — 1 
states are untrapped. 

In the present paper, we carry out scattering calcu- 
lations on Mg -f NH as a function of collision energy 
and magnetic field. The calculations use the MOLSCAT 
package [ITII , as modified to handle collisions in magnetic 
fields [18l|. The collisions that are of most interest are 
those of NH molecules that are initially in the magnet- 
ically trappable mj = +1 state, which may undergo in- 
elastic collisions to untrapped states with mj = and 
— 1. The calculations are carried out in a partly cou- 
pled basis set \nsjmj)\LMi,), where L is the end-over- 
end rotational angular momentum of the Mg atom and 
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the NH molecule about one another and Ml is its projec- 
tion on the axis defined by the magnetic field. Hyperfine 
structure is neglected. The matrix elements of the to- 
tal Hamiltonian in this basis are given in ref. p^ . The 
only good quantum numbers during the collision are the 
parity p = (— and the total projection quantum 
number M = rrij -f M^; the total Hamiltonian is block- 
diagonalized and scattering calculations are performed 
separately for each parity and M . The calculations in the 
present work use basis sets with rimax — 6 and Lmax — 8. 

MOLSCAT constructs a set of coupled equations in 
the intermolecular distance R and propagates them 
by the hybrid log-derivative method of Alexander and 
Manolopoulos [l^, which uses a fixed-step-size log- 
derivative propagator in the short-range region (2.5 < 
i? < 50 A) and a variable-step-size Airy propagator in the 
long-range region (50 < i? < 250 A). The log-derivative 
solutions for each M. are then matched to asymptotic 
boundary conditions [20| to obtain the scattering matrix 
S-^P and T-matrix T^p = / - S^p . It is useful to 
decompose the integral cross sections between NH lev- 
els (I a) = \nsjmj)) into sums of partial cross sections 
characterized by the L quantum number in the incident 
channel, 

L _ \ " \rpMp |2 

" MpL' 

where Ml = M - rrij, M'j^ — M - m'p ka is the 
wave vector for incoming channel a, with collision energy 
E = h'^kl/2ii, and fi is the reduced mass of the collid- 
ing system. Since we focus here on transitions among the 
n = 0, j — 1 levels, we abbreviate the labels a to just rrij . 
For molecules initially in the nij = +1 state, the most 
important quantities are the elastic and total inelastic 
cross sections; the latter is the sum of the state-to-state 
inelastic cross sections to the rUj = and nij = — 1 states. 

The mechanism of spin relaxation for ^Ti molecules has 
been studied extensively [111, [H, [H [13, [111. At suffi- 
ciently low energy, incoming channels with L > are sup- 
pressed by centrifugal barriers. The heights of the bar- 
riers are approximately E^f = {hL{L + l)//x)2 (54C6)2 , 
which is 23 mK for Mg-|-NH with L = 2. For L = (s- 
wave scattering). Ml — and hence M — rrij. Since M 
is conserved there is no outgoing channel with L' = for 
m'j = or — 1 and the dominant relaxation channels for 
s-wave scattering have L' = 2 in the outgoing channels. 

The coupling between channels with different mj oc- 
curs via the interplay of the spin-spin interaction and 
the potential anisotropy. The spin-spin term in the NH 
Hamiltonian mixes the n — and n = 2 states with the 
same j and rrij , and the potential anisotropy then mixes 
states of different L such that Anij + AMl = 0. 

In the absence of a magnetic field the thresholds for 
different values of mj are degenerate. The presence of 
centrifugal barriers in the outgoing channels strongly sup- 
presses the inelastic transitions, but the spin-relaxation 
cross section is nevertheless nonzero at finite energy 
[13, H^. Application of a magnetic field removes the 
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FIG. 1; (color online) Mg-|~NH elastic cr,^y£+i^+i (solid, red) 

and total inelastic cr,^j£:+i_^o + ""mj^+i^-i (dashed, blue) s- 
wave cross sections as a function of collision energy for various 
magnetic fields. Integral cross sections including p, d and f 
waves {L — 1, 2 and 3) are also shown for the elastic cross 
section (red, dots) and the total inelastic cross sections at 1 
and 10 G (blue, dots). This gives convergence for energies up 
to 100 mK. All cross sections are in A^. 
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FIG. 2: (color online) Mg-|-NH total inelastic s-wave cross 
sections a^Z+i^o + '^m^=+i^-i as a function magnetic field 
strength for various collision energies. All cross sections are 
in A2. 



degeneracy, increasing the kinetic energy in the outgoing 
channels and reducing the centrifugal suppression. By 
contrast, the elastic cross section is dominated by AL — 
collisions so are almost field-independent. 

Fig. [1] shows the s-wave elastic cross section ct^^^.^q 
and the total inelastic cross section cr^^.-^g for initial 
ruj ~ +1 as a function of energy for varying magnetic 
field strengths, while Fig. [5] shows the s-wave total in- 
elastic cross sections as a function of magnetic field for a 
number of different collision energies. It may be seen that 
the total inelastic cross section in the ultracold regime 
decreases dramatically as the magnetic field is reduced. 

For magnetic fields small enough that the NH Zeeman 



splitting does not exceed the centrifugal barrier height 
in the outgoing channel, the threshold behavior of the 
inelastic cross sections as a function of energy and mag- 
netic field can be understood as described by Volpi and 
Bohn [2lj. Applying the first-order distorted- wave Born 
approximation and approximating the scattering wave- 
function with the small-argument limit of spherical Bessel 
functions, they obtained a simple one-parameter formula, 

(2) 

where the factor cr^^/ is independent of energy and mag- 
netic field and AmjgiJ,oB is the linear Zeeman shift where 
Arrij = JTij — m'p g is the electron g-factor, and /ig is 
the Bohr magneton. For the centrifugal barrier to be ex- 
ceeded in the dominant outgoing channel (m^ — ~l, L' — 
2) for a zero-energy collision requires a magnetic field of 
86 G. In figure [T] it can be seen that the s-wave inelastic 
cross sections (i = ^ L' = 2) behave as 0. When 
the collision energy is less than the Zeeman shift, E^^2 
dominates and the s-wave cross section is proportional 
to E~2. For higher collision energies, the second term in 
(121) also contributes and the s-wave cross sections for low 
fields show E^ behavior. 

At low enough magnetic field, all the s-wave inelastic 
cross sections in figure [2] flatten out to a zero-field value 
proportional to E^. At higher field they enter a region 
of B^/^ dependence. This continues until the centrifugal 
barrier is exceeded in the outgoing channel and the cross 
sections then flatten off with a value proportional to E~ 2 . 

As mentioned in the Introduction, the general rule of 
thumb for sympathetic cooling to work is that the ratio 7 
of elastic to total inelastic cross sections must be greater 
than about 100. It can see in figure [T] that for small 
magnetic fields and low collision energies 7 is well in ex- 
cess of 100. However, at collision energies above ^ lO""* 
K, higher partial waves start contributing significantly to 
the total cross sections. The total cross sections incor- 
porating additional p, d and f partial waves (i = 1, 2 
and 3) are included in Fig. [T] for 1 and 10 G. There is a 
sharp peak in the d-wave inelastic cross section around 
75 mK, but everywhere else 7 remains in excess of 100 
until partial waves with L — A become important above 
100 mK. 

In order to assess the prospects of sympathetic cooling 
for NH with Mg, Fig. [3] shows a contour plot of 7 as a 
function of collision energy and magnetic field strength. 
At the top left of this Figure, inelastic collisions are too 
fast for sympathetic cooling to succeed. However, in 
an unbiased magnetic trap with zero field at the center, 
trapped molecules in rrij = -1-1 states at temperature T 
will be distributed according to a Boltzmann distribution 
with density p given by 

p/po=cxp(^^-^j. (3) 

At any given temperature on the energy axis of Fig. 
131 only about 0.1% of molecules will experience fields 
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FIG. 3: (color online) Contour plot of the ratio 7 of elastic 
to total inelastic cross section as a function of magnetic field 
strength and collision energy. The red line shows the max- 
imum field sampled by trapped molecules in the = 
state (B = QkT/gpo). 



greater than B = 6kBT/gfiQ, which is shown as a red 
line in Fig. [3l Most molecules in a magnetic trap will 
experience fields below this line. 

A precooled sample of molecules would initially have a 
temperature on the order of tens of milliKclvin. For tem- 
peratures up to about 10 mK, 7 is always greater than 
100 in the thermally allowed region. Even above this, 
7 remains close to 100. Nevertheless, it will be impor- 
tant to precool the molecules as much as possible before 
sympathetic cooling begins. As the sample is cooled to- 
wards sub- milliKclvin temperatures, the maximum mag- 
netic field strength sampled by the molecules decreases, 
7 increases and the trapped NH molecules become in- 
creasingly stable to collisional spin-relaxation. For an 
unbiased trap, 7 > 10^ at temperatures below 0.1 mK. 
However, the addition of a bias field to maintain the ruj 
quantum number would limit the increase in 7. 

Scattering at low energy depends strongly on the de- 
tails of the potential energy surface. The potential energy 
surface used in the present work is probably accurate to 
about 5%. To explore whether uncertainty in the poten- 
tial surface affects our conclusions, we consider the effect 
of a scaling factor A that produces a modified poten- 
tial energy surface I/'''='^'<=^ (i?, 9) = \V{R, 9). The s-wave 
elastic and total inelastic cross sections are shown as a 
function of A in Fig. [4] for a collision energy of 1 /xK at 
a field of 10 G. Both the elastic and inelastic cross sec- 
tions show strong resonance structures as Mg-NH bound 
and quasi-bound states cross the low-field-seeking thresh- 
old as the potential is varied. However, away from the 
strong resonant structures the ratio of elastic to inelastic 
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FIG. 4; (color online) S-wave elastic (solid, red) and total 
inelastic (dashed, blue) cross sections as a function of the 
potential scaling factor A, calculated for a collision energy of 
10~^K at a magnetic field of 10 Gauss. All cross sections are 
in A^. 



cross sections remains large. This indicates that the our 
conclusions are reasonably independent of the details of 
the potential energy surface and confirms that Mg is a 
good candidate for sympathetic cooling of magnetically 
trapped NH. 

In conclusion, we have carried out calculations on spin- 
changing collisions for NH('^I]^) molecules colliding with 
Mg atoms. We find that the ratio of elastic to inelastic 
cross sections exceeds 100, the factor required for sym- 
pathetic cooling to succeed, for a wide range of collision 
energies and magnetic fields. If precooled NH molecules 
at a temperature around 10 niK can be brought into con- 
tact with laser-cooled Mg, there is a good prospect that 
sympathetic cooling will succeed. The inelastic losses de- 
crease even further as the temperature decreases, so that 
once sympathetic cooling begins it will continue. 
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gramme EuroQUAM. 
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